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Transport through two quantum dots laterally embedded in Aharonov-Bohm interferometry with 
infinite intradot and arbitrary interdot Coulomb repulsion is analyzed in the weak coupling and 
Coulomb blockade regime. By employing the modified quantum rate equations and the slave- 
boson approach, we establish a general dc current formula at temperatures higher than the Kondo 
temperature for the case that the spin degenerate levels of two dots are close to each other. We 
examine two simple examples for identical dots - no doubly occupied states and no empty state. In 
the former, completely destructive coherent transport and phase locking appear at magnetic flux 
$ = $0/2 and $ = respectively; in the latter, partially coherent transport exhibits an oscillation 
with magnetic flux having a period of $o- 

PACS numbers: 73. 21. La Quantum dots, 73.23.-b Electronic transport in mesoscopic systems, 73.23.Hk 
Coulomb blockade and single-electron tunneling. 



I. INTRODUCTION 

Quantum dot (QD), a tiny engineered device accommodating a single electron or a few ones in three-dimensionally 
confined space, acts not only as a crucial ingredient for the realization of solid state quantum computation but 
also as a convenient tool to explore the effect of strong correlation manifested by discrete energy levels. Very rich 
phenomena, such as resonant tunneling, Coulomb blockade and the Kondo effect, arise in different circumstances 
to intrigue experimentists and theorists. To address the phase coherence of the transport, the QD or QDs are 
embedded in various Aharonov-Bohm (AB) geometries. Till now single do1ii*2iM£ or double dots&Z^ in two-terminal 
AB interferometer have been realized in experiments and the current oscillation of magnetic flux has been observed. 




FIG. 1: Schematic of a double-dot system in a parallel configuration between two leads in the presence of Aharonov-Bohm 
magnetic flux. 

In this paper, we consider two quantum dots parallely connected to two metallic reservoirs with magnetic flux 
penetrating the enclosed area. Due to the complexity of this system, most of the previous research concentrated on 
the noninteracting case to give the exact and general results of transport and to find effects of magnetic flux, position 
and difference of QDs energy levels, and band widths of coupling strength2ii2iiiii2ii&i!ii£. There were only a few 
thoretical studies of the interacting systems, which focused either on intradot correlation with spin for two identical 
dot a 12 i 16 i 17 , or on interdot correlation without spin for two different dot a 18 i 19 , or on the case having large intradot 
and interdot correlation without spin for two identical dotsi^. RefiS& is an exception where the intradot and interdot 
Coulomb repulsion, spin configuration and disparity of two dot levels were discussed at zero temperature covering the 
strong coupling Kondo regime. The present paper also deals with intradot and interdot Coulomb correlation, spin 
configuration and level disparity but is concerned with the Coulomb blockade regime at temperatures higher than the 
Kondo temperature. 

For the description of quantum transport through double-dot system, the "classical" rate equations must be modi- 
fied for nondiagonal density matrix elements responsible for transitions between isolated quantum states2ii22i2i22i2£. 
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Recently a modified quantum equations have been derived^ 6 - to study the transport of an interacting system using the 
slave-boson technique introduced by Zou and Anderson^ and incorporating the nonequilibrium Green's functions. 
The solutions are equivalent to the lowest-order gradient expansion, which is a good approximation for sequential 
resonant tunneling 28 ", and essentially accordant with previous analyses. Utilizing this kind of quantum rate equations, 
one is able to discuss the Coulomb correlation effect at arbitrary temperature. In the present paper, we apply this 
method to analyze a double-dot AB interferometer with interdot and intradot Coulomb repulsion, considering the 
lowest-order of the dot-lead coupling strength in transport. 



II. FORMULATION 



The Hamiltonian of two tunneling coupled quantum dots parallely connected to left and right leads with the presence 
of magnetic flux is described by the genetic tunneling model: 



H = egkaV^Cgka + ei Y. q^cia + Uni^mi + €2 c\ a C2a + Uni\nij + U' n\ a nz a i (1) 

other c cr aa' 

+ ^{tLloc\ ka C la + t Rla C Rka C lr7 + h.C.) + ^pL2ac\ ka C 2< y + t R 2aC Rka C 2 a + H.C.), 



ka ka 



where cL 2 , (ci(2) CT ) and c aka (c aka ) are creation (annihilation) operators of electrons in the dots 1(2) and in the left 
and right leads (a = L, R) with spin a. Each dot has a single spin degenerate orbital level £1(2) an d an infinite on-site 
Coulomb repulsion U and simultaneously there is an arbitrary finite interdot electrostatic correlation U' between 
them. We only consider the two dot levels are very close to each other, e 2 (i) = e,j±e/2, with a small variation e. 
The effect of AB flux $ is taken into account in the tunneling amplitude t al ( 2 )a by tp = 27r<I>/<I>o with the flux 
quantum $ = h/e. In the Peierls gauge, one generally chooses t* Lla — t^a — t* R2a = t R i a = |t|e' y/, 410 i 12 . The ath 
lead is supposed to be Fermi liquids in equilibrium state and has the same coupling strength function with two dots 

T aa {w) = 2TTj2ka |*al(2)<r| 2 <K W - £ a ka)- 

According to the slave-particle approach originated by Zou and Anderson, we introduce auxiliary operators 
e K fta [f2a)i d \a< to stand for the possible states of two dots as a whole: empty state |0) 1 |0) 2 , singly occupied 
state |cr) 1 |0) 2 (|0) 1 |cr) 2 ), and doubly occupied state |cr) 1 |(j') 2 respectively. In the slave-boson representation, the elec- 
tron operators of each dot are substituted by the slave-boson operators e^, and the pseudo-fermion operators fl a , 

ft 26.27 . 
J2a T1 ' ■ 

Cla = 6 1 fla + Y fL> d <™'> C 2<r = ^ho + Y fL'da'a, (2) 
a' a' 

with the completeness constraint e^e + ^ a .(/i ff /i £ r + f^ahv) + J2aa< d aa' d ^' = 1, for which these operators must 
be correctly quantized to make the sum rule for the physical electron valid and the commutators between them are 
satisfyin g 26 ; 29 : 

ee^ — 1, d aicr2 d cr , ia , 2 S aia ' 1 S cr2a i 2 , fiafj a ' = 3ij$aa' > (3) 

ed la' = e fL = /io-e f = U<r d a 'a" = d °°' e] = d °'a"fh = °- 

So the effective Hamiltonian is written in terms of these instrumental state operators: 

H eff = Y ^kJ aka c aka + ex Y flf la + e 2 ^ fLfz* + ( £ i + £ 2 + U') Y dl a ,d aa , (4) 

aka a a aa' 

+ Y} tL1 - c LkM ~h° + Y ft' d -') + tm^ Rka (^fia + Y fl' d °°') + h - c -\ 

ka a' a' 

+ Y^^LkA^ha + Y fL> d °'°) + t™J Rk Mh° + Y fl> d <r'") + h - c -i 

ka a' a' 

and the density matrix elements are expressed as poo = 1 0) x 1 0) 22 (0 1 1 (0 1 = e^e, pn a — |(t) 1 |0) 22 (0|i((t| = fl a fia, 

P22a = |0) 1 |£7) 22 (o-|l(0| = /^jW, Pddaa' = \ °) X W) 2 2 W 1 1 ( CT I = d la' d ™'' and Pl2a = |0) j |cr} 22 (0 1 1 ((7 | = f 2 Jl a , 

P2i CT = | ( r) 1 |0) 22 (a| 1 (0| = / 1 t CT / 2(7 . 
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Supposing the left and right leads are made from the identical material and the effective coupling strength is 
constant for the energy range of interest ^ ff (w) = T^^. Starting from the equations of motion of the density 
matrix elements and after the Langreth analytic continuation to decouple the dots and the leads interaction terms, 
the rate equations are acquired in the wide band limit: 

Poo = g J du, £L4 ff (G> lff + G> 22a ) + C a {G< la + G< 2a ) + B a G> X2a + B;G> 21r7 + D a G< 2a (5) 
Pn* = ^ / duj(A a G> la + l/2B a G> 12a + l/2B* a G> 21a + C a G< lla + l/2D a G< 2a + l/2D* a G< 21a ) (6) 

_ TZ J dLJ E [^'( G ll<xW + G lla'a"a) + B °' G d2la'a"a + B a lG d\2aa' a" + G °' ( G lla'aa>> 
J a'a" 

+ ^Lw) + B <?' G d21a'a"a + ^*a> G dl2aa'a" \ ' 

P22a = ^ j du{A a G> 22a + l/2B a G> 12a + l/2BlG> 2la + C a G< 22rT + l/2D a G< 12a + l/2D* a G< 2l(J ) (7) 

~ 4^ / diJ y^A A ^'( G 22a'aa" + G 22a'a"a) + B °' G d21a> aa" + B l' G dl2a"a'a + G a' ( G 22a> aa" 
a'a" 

+ ^22a'a"a) + At' ^d21a'aa" + At' ^dl2a" a' a\ ' 

Pi2« = ^ J M2A a G> 12a + Bl{G> lla + G> 22(7 ) + 2C a G< 2a + Dl{G< la + G< 2a )} (8) 

47T . 

+ G d21aa'a") + D a'( G dlla>a"a + G d22aa'a")] + *( £ 2 ~ ^l)Pl2<r, 



■ / (Lj [Aa'( G d21a'a"a + G d > 2\aa'a") + B a> [ G d\la> a" a + G d > 22aa'a") + G a' ( G d21a'a"a 

a'a" 



Pddaa - J dw^2[Aa{ G dllaaa> + G dllaa'a + G d22aaa> + G d22aa'a) + B ° ( G d21aaa> + G d21aa'a) ( 9 ) 



+ B a( G dl2a'aa + G d\2aaa' ) + G <?( G dllaaa' + G dllaa'a + G d22aaa' + G d22aa'a) + ^d G d2h 



aaa' 



PAfcff = ^ / ^E^^ll*-' + G rfH^' ff ) + 4x(G& 2<rff<r , + G^ Wff ) + B a G "d2\aaa' (™) 
a' 

+ B aG d21era , a + B a G dl2a'aa + B a G dl2aaa> + ^(^dHow' + G dllaa'a) + ^ ( G d22aa a' 
+ G d22aa'a) + B <? G d21aaa' + D^dZXaa'a + ^a G dl2a'aa + B a G dl2aaa' ] ' 

here A ff = / L (^)r LCT + Mu)^, C a = [1 - /t(w)]r L<T + [1 - /hM]^, B ff = / L (^)r LCT e-^/ 2 + / JJ (o;)r ilff e^/ 2 , 
£><7 = [1 - Sh^Wha^ 1 ' 2 + [1 - /i ? (w)]r jR(T e l¥ ' /2 , with the Fermi distribution function = l/^"-^ + 1), ^ 

the chemical potential of the ath lead at temperature 1//3. 

For a symmetric system T^o- = r^ CT = r CT , the dc current of the stationary state can be symmetrized I = {Il—Ir)/2 
and I a (ot = L, R) is obtained from the average of time derivative of electron number operator of the ath lead: 

1 = 1/ ^E^^H-^HH^n. + G^ (ii) 

a a'a" 
- G dllaa'a" ~ G d22aa'a")\ + C0S 77 L/i M ~ f R^)] [ G el2a + G e21a - G el2a - G e21a + E ( G d21aa'a" 

+ G> dl2a«aa> G d21aa>a» G dl2a»aa>)\ +ism^[f L (uj) + f R (uj)} [G> 2 1<7 G> 2a + ]T {G> 12a „ a<J , 

a' a" 

- G 'd2.\aa' a" )] + 2* tin |[1 - M^ + Mf)] [G < 1(T - G< 12(7 + ^ (G< 2ff , w - G^< Wff „)]}- 
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This expression is similar with that of Ref^ where the interdot correlation was not included. The Green's functions 
quoted in the above are defined as follows: 

G e ua = « e\t)f ia (t)\fl(t')e(t') », Geijcr =« e* (t)f icr (t) \f} a (t')e(t') », (12) 
G d XXa>aa» = « fl(t)d^ (t) |4„ a , (t')fxa" (*') », G^ lw , a =« fl„ (t)<W (*) |4 v (i')/la (*') », 
G d 22a>aa» = «fL(t)d r r',(t)\dl <J „(t')f 2a „(t') »,G d22aa , a „ =« fl(t)d aa , (t) |4, CT „ (f)f 2 a» (*') », 

G txixcr' (j" g = «ft„(t)d^(t)\dl a (f)f 1(T (t') »,G' d21aa , (T „ =«fl{t)d aa ,{t)\dl l(rl {t')h all {t') », 



G. 



dlla'oo" 



= « fLW«>«{t)\dl,,At')h°"{t') »,Gdi2**'*» =« fl(t)d^(t)\dl a „(t')f 2a „(t') » 



In weak coupling approximation and with small level discrepancy e, the correlation Green's functions in the isolated 
two dot system are: 



G<»A" 

<0 
eija 



G dXXa'ac" ( W 
Gd^a'crcr"^ 



'<0 
dXXa 

'<0 

d22aa' a 
-i<0 



'a-a^ 



G'd2Xa>a"a( U} 



G 



'<o 

d21crcr' cr 



•A" 



GdTLa'crcr" ( U 

r<<o 



>A»> 



2mp iia 8{u) - e d ),G< lcr (ui) = -2irip o8(uj - e d ), 
2nip ij(T 8(u) - e d ) 7 G>° a {uj) = 0, 

5 aa »2mpddaa'5{uJ -e d - U'), G^ lrT , a(7 „(uj) = -8 aa »2mpu< 7 8{u -e d - V), 

d <T<7 »2'Kipdda'aS(0J -e d ~ U'), G^Lw' ( W ) = -$aa"2irip22aS(u ~ Cd ~ U'), 

Saa'8<r<T"2Kipdda<T8(u - e d - U'), G d \® a , a „ a {oj) = -8 aa >8 aa »2mp lla 8(uj -e d - V), 
8 aa '8 aa n2%ip d d <TtT 8{w -e d - V), Gj$b , „(w) = -8 aa >8 aa »2Trip22a8(uj -e d - U'), 



(13) 



0,G>° 



d22ao" 

^d2Xa'a"a( U ) = ~8aa"^ipX2aK^ - € d ~ U'), 

0, G>l ala „{Lo) = -8 rT<T n2mp 1 2 a 8(u; - e d - U'), 
0,G^°„,„„„(u) = -8 aa '8 aa »2mp 12 a&{u W), 



'i u d2Xa' 'aa" 

n r ,>0 



(w) = -<5 - -'(5 - -''27rip2i -^(w - ed - U'). 



Note that by using the lowest-order gradient expansion with slowly varying in the center-of-mass time T and rapidly 
varying in the relative time t and after the Fourier transformation from t to u>, the same results as the above can be 
acquired. 

Inserting Eq.(13) into Eq.(5) — (10), we get the final quantum equations: 

POO = 2q W00 + beta\a(j)\i a + p 2 2a) + 02aPX2a + PtaP^Xa], ( 14 ) 

a 

PXXa = "lcr/OoO - {$Xa + &Xa>)pXXa - l/2(/?2 CT + a 2a )px2a - l/2(/3| CT + a 2a )p2Xa + flla'Pddaa', 

a' cr' 

p22a = &laP00 - (fl\o + ^ "icr' )/°22cr - l/2((3 2 a + a 2a )PX2a ~ V 2 ($2<7 + a* 2a ) P2Xa + fixa'Pdda'a, 

cr' cr' 

P21o- = "2(7/000 - l/2(/3 2ff + a 2o -)(pil C r + p22<7 ) - (/3lcr + fii^ ) ,02 lcr + fcaPddaa + «(ei ~ e 2 )p21 CT , 

Pcidcro- = Q!1(t(P11(T + P22cr) + &2aPX2cr + &2 a P21cr ~ 2(3 lcr Pddaa , 
PddaS — aiaP225 + OiXaPXXa ~ {Pier + Pxa)pddaa- 



Here. 



"lcr = fL(ed)T La + fR(e d )T R<T ,f3 la = [1 - /i(e rf )]r iCT + [1 - /K(ed)]r flcr , (15) 

aia = h{e d + U')T L<y + f R (e d + U')T Ra J ll7 = [1 - f L {e d + U')]T La + [1 - f R (e d + U')]T Ra , 

a 2 a = f L (ed)r La e-^/ 2 + f R (e d )r Ra e iv/2 , ha = [1 - h(e d )]I ' Ltr e~^ + [1 - f R (e d )}T Ra e iv/2 , 

&2a = fL{e d + U')T La e-^' 2 + f R (e d + U')T Ra e^' 2 , 

ha = [I - fL(e d + U')]Y La e- 1 */ 2 + f R (e d + U')]T Ra e^l 2 . 



Complemented with the completeness relation p o + 2p\\ + 2p 2 2 + 2p dd + 2p dd 
be solved to determine the expectation values of the matrix in steady states. 



1, the closed equations Eq.(14) can 
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III. DISCUSSION 



In the equilibrium state, only the diagonal distribution probabilities are nonzero, poo — l/Z, p\\ a = 



-Pea 



I Z, Pddo 



_ (S -/3(2e d + C/') 



/Z,p 2 2a = 



/Z and Z = 1 + 4e /3cd + 4e / 3 ( 2e <i+ c/ ) f w hich can be readily gotten from the above 



equations set with f L (e d ) = f R {e d ) = (e* 3 ^ + l)-\f L {e d + W) = f R {e d + U') = { e ^+ u '^ + l}" 1 supposing the 
equilibrium chemical potential pl = PR = 0. The results meet the classical Boltzmann distribution in weak coupling 
limit. It is obvious that the occupation number of dot 1 with certain spin c\ a c\ a — p\\ a + p dd aa' has nothing to 
do with the magnetic flux $ and the effect of dot level disparity e is not taken into account, which is the same as dot 
2, and in the spin symmetry space p ddaa = p ddaS . 

Applying external voltage on the two leads, the whole system is driven out of equilibrium and off-diagonal elements 
play a vital role in transport, whereupon some novel features arise. In the spin symmetry space a{a)n2) a , P{P)i(2)a, 
r l(r)<j are all independent of spin, so for simplicity we subtract their spin index a. After straightforward derivation, 
the solutions for these closed equations and the final formula of current in steady states are obtained in Appendix. In 
the following we will see that the occupation number of the dot 1(2) shows explicit AB oscillations originating from 
the interference of two dot states p\2a and P2ia, which is certificated by the existence of phase related coefficients 
«2(t, fca, &2tn P20 in the closed equations only stemming from the existence of p\2a and p2\a- Furthermore, p dda a is 
distinct from p dda g. This difference can be understood in this way. Interference between electrons in the leads and in 
the dots survives only when both dots are empty, or one dot is empty and the other is occupied by or obth are occupied 
by the same spin as the electron coming from the lead. So the bias-induced coherence of two dot states p\2a and p2\a 
are only associated with spin diagonal double-occupied state probability pddaa and not with spin nondiagonal one 
Pddaa, otherwise the path the incoming electron travels can be identified. This is distinct from double-dot in scries 
where the path is a definite one including both dots 2 ^. From the equation of pddaa, we can see that effect of p\2a and 
P2\a is to accelerate its change with the time. 
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FIG. 2: Left: AB oscillations of dot population «i(2) with interdot correlation U' — and energy disparity e — at temperature 
ksT = 20 and bias voltage eV = 50 for different average level e^: —l(solid), — 1.5(dash), — 2(dot), — 2.5(dash dot). Right: AB 
oscillations of current for different energy disparity e with interdot correlation U' = 20 and average level e<j = — 1 at temperature 
k B T = 20 and bias voltage eV = 50. e = H = F = 1. 



Note that the temperature is assumed higher than the Kondo temperature so that the Kondo correlation which is 
attributed to the strong coupling between leads and dots can be negligible. In the sequential tunneling picture, the 
tunneling contribution to transport is up to first order in T in the Coulomb blockade regime. We can directly calculate 
the current from Eq.(2I) which is valid for arbitrary bias voltage eV , interdot repulsion [/', and dot energy position 
e d in broad ranges and for tiny level variation e between two dots. For convenience, the symmetric bias voltage 
PL = —PR. = eV/2 is applied, then all the results are invariant under the reversals of magnetic flux tp— > — ip and of 
bias voltage eV^ — eV simultaneously. This is in accordance with the Onsager relation for two-terminal setups^. In 
general, the occupation number of dot 1(2) and the current through the system all have AB oscillations, combinations 
of cos ip and sine/?, with a period of $0, and the amplitudes depend on the resonant energy level e d and e d + U' 
relative to the chemical potentials of two leads and small disparity amplitude e. It is also supposed that ksT^>T, 
I ei I , I £2 1 , and r^>|e|, so the lowest-order transport dominates. In Fig. 2 we plot AB oscillations of the population in 
dot with different level positions and the current with different energy discrepancys respectively. For two identical 
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FIG. 3: Current versus interdot correlation U' with various temperatures for average dot level = —1, energy disparity e = 0, 
and bias voltage eV — 50 at tp = (left) and at tp = n (right), e = h = T — 1. 



dots e = 0, all oscillations are symmetric about tp = including only cos if with larger amplitudes for deeper levels. 
The population reaches maximum at tp = and minimum at if = ir. For two distinct dots e^O, the symmetry about 
ip = is lost, which gets more apparent with larger e for oscillations following simp take effect, and the extremes at 
(p = and ip — 7r also disappear. Fig. 3 shows current versus interdot correlation with various temperatures for tp = 
and tp = 7T. We can see that at finite temperature and finite bias the interdot correlation does not always inhibit 
transport. Only when U' is large enough, the prohibition effect appears, especially at tp = it, and this is mitigated by 
raising the temperature. 

In the following two specially simple cases with equal energy levels of two dots are under consideration. First 
is that doubly occupied states are forbidden by infinite interdot Coulomb repulsion U', such that pddaa' — 0, and 
Q'i = 0, a-2 = 0. The steady solutions are 

(2 - f_L - f R )[(l - f L )fn + (1 - fn)f L ] nfil 
Pl1 2(f L + f R -3f L f R +l)(2-f L -f R ) ' 1 } 

fh - Ir r fit - Jl <P . . . <Pi 

Pl2 = 7777 7 777 — 7 7T Y7, 7 7~ cos 77 + 1 sm " 



POO = 



2(/ £ + f R - 3/i/fl + 1) L 2 - f L - f R 2 2 
Jx/g ~ /l — /g + 1 

/L+/fl-3/i/ fl + l' 



and the current is 

r 27rer 2(f L -f R )(l-f R )(l-f L ) 



h (h + f R - 3/i/fl + 1)(2 - f L - f R ) 



(1 + cos^). (17) 



It is evident that in this situation transport through the whole system is totally coherent, even though the population 
of dot 1(2) shows no AB oscillations. With tp = 0, ±27r, ±4-7t, the current reaches the maximum, and phase locking 
always happens for any voltage and dot levels. With tp — ±7r, ±37r, ±57r, the current vanishes and the interference 
between two dots is completely destructive, since the whole system now is equivalent to the sum of two separated 
subsystems and no transport will take place for any bias, which conforms to Ref^. Considering only the first order 
of r in sequential tunneling regime, we get the full coherence. 

Second is that empty state is forbidden in deep level status with far below the Fermi level and td + V just above 
it in equilibrium state, so that poo = 0, and j3\ = 0, /?2 = 0. The steady solutions are 

(2 - ft - fn)[(l - hWfii + h~ fneostp) + (1 - f R )(2f L + f R - f L costp)} 

Pll — = = = = = , ( ±o I 

4(/l + / fl )(3 - cos^)(2 -f L - f R ) + 2(f L - f R f{l - costp) 

P12 = Af ~ l ~ . [(2 f l fa) cos | + i sin | (/ l - f R ) - (4 cos g - 2z sin | - **f h ) 
4(/l + Jr) 2 2 2 2 2 - f L - f R 

(2 - /l - fn)[(l ~ }l){Vr + h~ /flcos^) + (1 - fn)(2f L + f R ~ hcostp)} 
4(/l + /a)(3 - cos(p)(2 - f L - f R ) + 2{f L - f R ) 2 (l - costp) 
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1 (4 - f L - f R ) [(1 - h)2f R + f L - f R cos tp) + (1 - /r) (2f L + f R - f L cos ip)] 

Pdd ~ 



Pdd 



2 4(/ L + /«) (3 - cos (2 - / L - + 2(/ L - f R ) 2 (1 - cos 

(/l + /«)[(! - !l){Vr + h - Ircos^) + (1 - f R )(2f L + f R - f L cos^)] 



4(/l + / fl )(3 - cos<^)(2 - /t - f R ) + 2(f L - / fl ) 2 (l - cosip) 
and the current is 

/ = ^ {1 + 2(f L + f R ) + (l- f L - f R ) cos <p - [(4 - / L - (19) 

,7 ,f x ^ 1 (1 - /l)(2/a + /l - /flcos^) + (1 - / fl )(2/ L + f R - / L cosy) 1 
+ (4 — S fL — o t R ) cos w\ = = = = = = k 

2(/L + /fl)(3-C0S V )(2-/ i -/ fl ) + (/ L -/ fl )2(l-C0S V ) 

which is partially coherent, and the AB oscillations for the population in dots and the current all follow the cos ip 
form. 



IV. SUMMARY 

To conclude, we have studied the quantum transport through a parallel double-dot structure with the Aharonov- 
Bohm magnetic flux infiltering the two path closed region. Infinite intradot and arbitrary interdot Coulomb repulsion 
have been considered by employing the slave-boson technique introduced by Zou and Anderson. In weak coupling 
and sequential tunneling regime we use the " classical" quantum rate equations combined with noncquilibrium Green's 
functions to determine the density matrix and to calculate the current flowing through the system of stationary state. 
The results show that external bias voltage induced superposition of two dots states is the cause of phase coherence of 
AB oscillations of population of each dot and the current as a function of magnetic flux having a period of $o- Two 
simplest cases are discussed as examples for two identical dots. Wc find that if there arc no doubly occupied states, 
the current is totally coherent - completely destructive interference exists with ip odd multiples of n and maximum 
always reaches with ip even multiples of tt; if there is no empty state, the population and the current are partially 
coherent and amplitudes depend on the resonant energy levels relative to the chemical potentials of two reservoirs. 

This work was supported by the National Science Foundation of China, the Special Funds for Major State Basic 
Research Project, and the Shanghai Municipal Commission of Science and Technology. 

V. APPENDIX 

The solutions of density matrix elements in stationary states are: p\\ a — p\\ s — pn, p 22a = P225 = P22, P2\a = 

P21S = P21, Pl2a = P\2S = Pl2, Pddcrcr = PddSS = Pdd, PddaS = PddSa = Pdd, 

Poo = (C 2 A 1 -C 1 A 2 )/(B 1 A 2 -B 2 A 1 ), (20) 

P11 = {CiBi-dBiy^AtBi-AiBi)], 

P21 = {{2a 2 - (3 2 ) Poo - [& 2 +/3 2 + {2 + &i/jh)ih]{.Pxi + P22) + &}/[2(/?i + 25i + ie% 

Pdd = [1 -poo - (2 + ai//3i)(pn +p22)]/2, 

Pdd = 61O11 +P22)/(2/3i), 

and p 22 = pn, p 12 = p* 21 . Here, 

A\ = 4(A 2 + e 2 ) - A(/ L + f R cos tp) [1 + f L - f L + 7(1 - f L )] - A(f R + f L cos <p) [1 + f R - f R 

+ 7(1 - /«)] + c sin - / L )(l + 7) + hfn - Ml], 
Bi = A[2/ L (/ L + /flcos^) + 2f R (f R + / L cos<^) - (1 - f L )(h + f R cos<p) - (1 - f R )(f R 

+ f L cos if)} + (2 - f L - f R ) (A 2 + e 2 ) - e sin <p[2(f R f L - f L f R ) + f L - f R ] , 
Ci = -(2 - f L - f R )(A 2 + e 2 ) + A[(l - f L ){h + f R cos<p) + (1 - f R )(f R + f L costp)} 

+ esinip(f L - f R ), 

M = (2 - f L - /k)(A 2 + e 2 ) - A(l - f L )[l + f L - f L + cos^(l + f R - f R ) + 7 (1 - f L + cos<^ 



8 



- COS lff R )] - A(l - f R )[l + f R - f R + COStp(l +f L ~ /l) +7(1 - fR + COSV3 - cosipf L )] 

+ esin^[(l - f R )(f L + 7 - jf L ) ~ (1 - + 7 - 7/*)]. 

£ 2 = A(2/ L + /i - 1)(1 - f L + cos <p - cos p/ R ) + A(2/r + - 1)(1 - /« + cos 99 - cos <^/ L ) 

- 2(/ L + / fi )(A 2 + e 2 ) - esm^[(l - f R )(l + f L ) - (1 - /z)(l + /«)], 

C 2 = A[(l - / L )(l - / £ + costp - cos(pf R ) + (1 - /h)(1 - f R + cos<p- cos <pf L )] 

- esin^[(l - f R )(l - ~f L ) - (1 - /l)(1 - /*)], 

and the current formula becomes 
2neT 

I = — [(h - / fi )(2poo + Pu + P22) + cos^/2(/ L +f L -f R - f R )( Pl2 + P21 ) (21) 
- i sin ^/2(/z, + /x + f R + f R - 2)(pia - P21) + 2(/i - / B )(pn + p 2 2 + Pdd + 
where A=2 - f L - f R + 2f L + 2f R , 7 =(4 - / £ - /«)/(2 - f L - f R ), 

and f HR) =f L{R) {e d ) = + j L{R) =f L{R) {e d + U>) = {e^^+^'-^c*)) + 1}-*. 
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